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We translate the action of local Clifford operations on graph states into transformations on their 
associated graphs - i.e. we provide transformation rules, stated in purely graph theoretical terms, 
which completely characterize the evolution of graph states under local Clifford operations. As we 
will show, there is essentially one basic rule, successive application of which generates the orbit of 
any graph state under local unitary operations within the Clifford group. 
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I. INTRODUCTION 



Stabilizer states and, more particularly, graph states, 
and (local) unitary operations in the Clifford group have 
been studied extensively and play an important role in 
numerous applications in quantum information theory 
and quantum computing. A stabilizer state is a multi- 
qubit pure state which is the unique simultaneous eigen- 
vector of a complete set of commuting observables in the 
Pauli group, the latter consisting of all tensor products of 
Pauli matrices and the identity (with an additional phase 
factor). Graph states are special cases of stabilizer states, 
for which the defining set of commuting Pauli operators 
can be constructed on the basis of a mathematical graph. 
The Clifford group consists of all unitary operators which 
map the Pauli group to itself under conjugation. As 
the closed framework of stabilizer theory plus the Clif- 
ford group turns out to have a relatively simple mathe- 
matical description while maintaining a sufficiently rich 
structure, it has been employed in various fields of quan- 
tum information theory and quantum computing: in the 
theory of quantum error-correcting codes, the stabilizer 
formalism is used to construct so-called stabilizer codes 
which protect quantum systems from decoherence effects 
Pj ; graph states have been used in multipartite purifica- 
tion schemes and a measurement-based computational 
model has been designed which uses a particular graph 
state, namely the cluster state, as a universal resource for 
quantum computation - the one-way quantum computer 
[3|; (a quotient group of) the Clifford group has been 
used to construct performant mixed-state entanglement 
distillation protocols ; most recently, graph states were 
considered in the context of multiparticle entanglement: 
in [3 the entanglement in graph states was quantified 
and characterized in terms of the Schmidt measure. 

The goal of this paper is to translate the action of local 
Clifford operations on graph states into transformations 
on their associated graphs - that is, to derive transforma- 
tion rules, stated in purely graph theoretical terms, which 



completely characterize the evolution of graph states un- 
der local Clifford operations. The main reason for this 
research is to provide a tool for studying the local unitary 
(LU) equivalence classes of stabilizer states or, equiva- 
lently, of graph states 0| - since the quantification of 
multi-partite pure-state entanglement is far from being 
understood and a treatise of the subject in its whole is ex- 
tremely complex, it is appropriate to restrict oneself to a 
more easily manageable yet nevertheless interesting sub- 
class of physical states, as are the stabilizer states. The 
ultimate goal of this research is to characterize the LU- 
equivalence classes of stabilizer states, by finding suitable 
representatives within each equivalence class and/or con- 
structing a complete and minimal set of local invariants 
which separate the stabilizer state orbits under the ac- 
tion of local unitaries. We believe that the result in this 
paper is a first significant step in this direction. 

In section IV, we will show that the orbit of any graph 
state under local unitary operations within the Clifford 
group is generated by repeated application of essentially 
one basic graph transformation rule. The main tool for 
proving this result will be the representation of the stabi- 
lizer formalism and the (local) Clifford group in terms of 
linear algebra over GF{2), where n-qubit stabilizer states 
are represented as n-dimensional linear subspaces of 
which are self-orthogonal with respect to a symplectic in- 
ner product and where Clifford operations are the 
symplectic transformations of Z2" |j, 0] ■ 

This paper is organized as follows: in section II, we 
start by recalling the notions of stabilizer states, graph 
states and the (local) Clifford group and the translation 
of these concepts into the binary framework. In section 
III, we then show (constructively) that each stabilizer 
state is equivalent to a graph state under local Clifford 
operations, thereby rederiving a result of Schlingemann 
jsj. Continuing within the class of graph states, in section 
IV we introduce our elementary graph theoretical rules 
which correspond to local Clifford operations and prove 
that these operations generate the orbit of any graph 
state under local Clifford operations. 
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II. PRELIMINARIES 



A. Stabilizer states, graph states and the (local) 
Clifford group 



Let Gn denote the Pauli group on n qubits, consisting 
of all 4 • 4" n-fold tensor products of the form a wi ® W2 ^ 
• ■ • (g) u„, where a £ {±1, ±i} is an overall phase factor 
and the 2 x 2-niatrices Vi {i = l,...,n) are either the 
identity ctq or one of the Pauli matrices 
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The Clifford group C„ is the normalizer of Gn in U{2"), 
i.e. it is the group of unitary operators U satisfying 
UGnU^ = Gn- We shall be concerned with the local Clif- 
ford group C^, which is the subgroup of C„ consisting of 
all n-fold tensor products of elements in Ci. 

An n-qubit stabilizer state \'ip) is defined as a simul- 
taneous eigenvector with eigenvalue 1 of n commuting 
and independent f3| Pauli group elements Mi. The 
n eigenvalue equations Mi\ip) = \'ip) define the state 
IV') completely (up to an arbitrary phase). The set 
S := {M € Gn\M\ip) = \'4))} is called the stabilizer of 
the state It is a group of 2" commuting Pauli op- 
erators, all of which have a real overall phase ±1 and 
the n operators Mi are called generators of <S, as each 
M e <S can be written as M = Mf ' . . -M^, for some 
Xi £ {0, 1}. The so-called graph states 2; 3j constitute 
an important subclass of the stabilizer states. A graph 
is a pair G — {V, E) of sets, where V is a finite subset 
of N and the elements of E are 2-element subsets of V . 
The elements of V are called the vertices of the graph 
G and the elements of E are its edges. Usually, a graph 
is pictured by drawing a (labelled) dot for each vertex 
and joining two dots i and j by a line if the correspondig 
pair of vertices G E. For a graph with \V\ — n 

vertices, the adjacency matrix 9 is the symmetric binary 
n X n-matrix where 9ij = 1 if G E and 9ij ~ oth- 

erwise. Note that there is a one-to-one correspondence 
between a graph and its adjacency matrix. Now, given 
an n- vertex graph G with adjacency matrix 9 one defines 
n commuting Pauli operators 



n „ 



where cri^ , a^y' , crt^ are the Pauli operators which have 
resp. cTj; , tjj, , tJ^ on the ith position in the tensor 
product and the identity elsewhere. The graph state 
|V'/iiM2 - p,i (G*)), where [li S {0,1}, is then the stabilizer 
state defined by the equations 



.(G)). 



Since one can easily show that the 2" eigenstates 
|''/'A'iA»2- .Mn (G)) are equal up to local unitaries in the Clif- 
ford group, it suffices for our purposes to choose one of 



them as a representative of all graph states associated 
with G. Following the literature j^, we denote this rep- 
resentative by \G) := |V'oo...o(G)). Furthermore, if the 
adjacency matrices of two graphs G and G' differ only 
in their diagonal elements, the states \G) and |G") are 
equal up to a local Clifford operation, which allows for 
an a-priori reduction of the set of graphs which needs 
to be considered in the problem of local unitary equiv- 
alence. The most natural choice is to consider the class 
O C of adjacency matrices which have zeros on the 

diagonal. These correspond to so-called simple graphs, 
which have no edges of the form or, equivalently, 

none of the points is connected to itself with a line. From 
this point on, we will only consider graph states which 
are associated with simple graphs. 



B. The binary picture 

It is well-known [3, 0, 0| that the stabilizer formal- 
ism can be translated into a binary framework, which 
essentially exploits the homomorphism between Gi, ■ and 
Z2,+ which maps ctq = <^oo ^ 00, ax = (Joi ^ 01, 
az — <Jio 1-^ 10 and ay — an i-^ 11. In Zj addition is to 
be performed modulo 2. The generalization to n qubits 
is defined by au^y^ (g) ■ ■ ■ (g) cr„„„„ = cr(„j...„^|.„j...^_^) i-^ 
(ui . . .Un\vi . . . Vn) S ^i" ; where Ui, Vi £ {0, 1}. Thus, an 
n-fold tensor product of Pauli matrices is represented as 
a 2n-dimensional binary vector. Note that with this en- 
coding one loses the information about the overall phases 
of Pauli operators. For now, we will altogether disregard 
these phases and we will come back to this issue later in 
this paper. 

In the binary language, two Pauli operators Oa and 
ab^ where a,6 G Z^", commute iff a^Pb = 0, where the 

2n X 2n-matrix P = ^ q defines a symplectic inner 

product on the space Z^". The stabilizer of a stabilizer 
state then corresponds to an n-dimensional linear sub- 
space of Z2" which is its own orthogonal complement 
with respect to this symplectic inner product. Given a 
set of generators of the stabilizer, we assemble their bi- 
nary representations as the columns of a full rank 2n x n- 
matrix S, which satisfies S^PS = from the symplec- 
tic self-orthogonality property. The entire stabilizer sub- 
space consists of all linear combinations of the columns 
of S, i.e. of all elements Sx, where a; G Zj. The ma- 
trix 5, which is referred to as a generator matrix for 
the stabilizer, is of course non-unique. A change of gen- 
erators amounts to multiplying S to the right with an 
invertible n x n-matrix, which performs a basis change 
in the binary subspace. Note that a graph state which 
corresponds to a graph with adjacency matrix 9, has a 



generator matrix S — 



Finally, it can be shown 



0, that, as we disregard overall phases, Clifford op- 



erations are just the symplectic transformations of 
which preserve the symplectic inner product, i.e. they 
are the 2n x 2n-matrices Q which satisfy Q^PQ = P. 
As local Clifford operations act on each qubit separately, 

\ A B 

they have the additional block structure Q = (j 

where the n x n-blocks A, B, C, D are diagonal. In this 
case, the symplectic property of Q is equivalent to stat- 

Aii Bi- 

ith qubit, is invertible. The group of all such Q will be 
denoted by C'. 

Thus, in the binary stabilizer framework, two stabilizer 
states and |?/'') with generator matrices S and S' are 
equivalent under the local Clifford group iff there is 
a, Q £ C'' and an invertible R £2,^ 



ing that each submatrix 



which acts on the 



such that 



QSR = S'. 



(1) 



Note that the physical operation which transforms 
into \i]y) is entirely determined by Q; the right matrix 
multiplication with R is just a basis change within the 
stabilizer of the target state. 



III. REDUCTION TO GRAPH STATES 



In this section we show that, under the transforma- 
tions S — > QSR, each stabilizer generator matrix S can 
be brought into a (nonunique) standard form which cor- 
responds to the generator matrix of a graph state. 

Theorem 1: Each stabilizer state is equivalent to a graph 
state under local Clifford operations. 

Proof: Consider an arbitrary stabilizer with generator 

' The result is obtained by proving the 



matrix S 



X 



existence of a local Clifford operation Q E C' such that 
Z' 

' has an invertible lower block X'. Then 



QS = 



X' 



S' := QSX'-^ = 



where Z'X' ^ is symmetric from the property S'^ PS' = 
0; furthermore, the diagonal entries of Z'X'^^ can be put 

11" 



to zero by additionally applying the operation 



1 



to 



the appropriate qubits, since this operation flips the zth 
diagonal entry of Z'X'^^ when applied on the zth qubit. 
Eventually we end up with a graph state generator matrix 
of the desired standard form. 

We now construct a local Clifford operation Q that 
yields an invertible lower block X'. We start by per- 
forming a basis change in the original stabilizer in order 



to bring S in the form 



S 



Rz Sz 
Rx 



such that Rx is a full rank n x /c-matrix, where k = 
rank X\ the blocks i?z, Sz have dimensions n x k, resp. 
n X {n — k). The symplectic self-orthogonality of the 
stabilizer implies that S^ Rx = 0. Furthermore, since 
Sz has full rank, it follows that the column space of Sz 
and the column space of Rx are each other's orthogonal 
complement. 

Now, as Rx has rank k, it has an invertible k x k- 
submatrix. Without loss of generality, we assume that 
the matrix consisting of the first k rows of Rx is invert- 
' R^ ' 

' where the upper k x /c-block Rl. 



ible, 



Rx 



Rl 

is invertible and R^ 
titioning Sz similarly in a k x 

{n — k) X [n — fc)-block Si , i.e 



has dimensions [n — k) x k. Par- 
{n — fc)-block Si and a 
S^ 1 

the prop- 



erty SJRx = then implies that S^ is also invertible: for, 
suppose that there exist x G Z,2~'' such that {S^^x = 0; 
then the n-dimensional vector v :— (0, . . . , 0, x) satisfies 
Sjv = and therefore v = RxU for some y G Z^. This 
last equation reads 



Ri 



Rb 
Rly 



Since R^. is by construction invertible, = implies 
that ?/ = 0, yielding x = Rl.y = 0. This proves the 
invertibility of S^.. 

In a final step, we perform a Hadamard transformation 

^ Q on the qubits fc-f 1, . . . , n. It is now easy to verify 

that this operation indeed yields an invertible lower nxn- 
block in the new generator matrix, thereby proving the 
result. □ 

This proposition is a special case of a result by 
Schlingemann |^, who showed, in a more general context 
of d-level systems rather then qubits, that each stabilizer 
code is equivalent to a graph code. 

Remark: overall phases - Theorem 1 implies that our 
disregard of the overall phases of the stabilizer elements 
is justified. Indeed, this result states that each stabilizer 
state is equivalent to some graph state \'4'^^l^l2■■■fJ.ni^))l 
for some /z^. As such a state is equivalent to the state 
|G), there is no need to keep track of the phases. 

Theorem 1 shows that we can restrict our attention to 
graph states when studying the local equivalence of sta- 
bilizer states. Note that in general the image of a graph 

I A B 

state under a local Clifford operation Q ~ CD 
not again yield another graph state, as this transforma- 



need 



tion maps 



'9' 




■ e ' 








I 


I 









(2) {e^') 



for g 9. The image in |5J| is the generator matrix of a 
graph state if and only if (a) the matrix C9 + _D is non- 
singular and (&) the matrix 9' := {AO + B) {C9 + Dy^ 
has zero diagonal. Then 



Q 



{C9 + D) 



is the generator matrix for a graph state with adjacency 
matrix 9' E Q. Note that we need not impose the con- 
straint that 9' be symmetric, since this is automatically 



the case, as 



is the image of a stabilizer generator 



matrix under a Clifford operation, and thus 



[e'^ i]p 



0. 



These considerations lead us to introduce, for each Q S 
C', a domain of definition dom((3), which is the set con- 
sisting of all G O which satisfy the conditions (a) and 
(6) . Seen as a transformation of the space Q of all graph 
state adjacency matrices, Q then maps 9 G dom((5) to 



Q{9) {A9 + B) {C9 + oy 



(3) 



In this setting, it is of course a natural question to ask 
how the operations Q affect the topology of the graph 
associated with 9. We tackle this problem in the next 
section. 

To conclude this section we state and prove a lemma 
which we will need later on in the paper. 



Lemma 

s.t C9 -I 
A 
C 

such that 



1: Let 9^0 and C,D be diagonal matrices 
- D is invertible. Then there exists a unique 

jj G C', where A, B are diagonal matrices, 

6 G dom{Q). 

Proof: Note that, since C9 -I- is invertible, we only 
need to look for a Q s.t. Q{9) has zero diagonal in order 
for 9 to be in the domain of Q. First we will prove the 
uniqueness of A and B: suppose there exist two pairs of 
diagonal matrices A, B and A' , B' s.t. 



G C' 



and 9 G dom(Q), 9 G dom{Q'). Denoting 9^ := A9 + B, 
9'^ := A' 9 + B' and 9^ := C9 -D, we have Q{9) = 9^9-^ 
and Q'{9) ^ 9',9-^. Now, denoting by zf,zf,xf the 



" A 


B ' 


M' ■■= 


'A' 


B'' 


C 


D 


C 


D 



rows of resp. 9z,9'^,9x, the crucial observation is that ei- 
ther zf = zf or zf = zf+xf for alH = 1, . . . , n, which is 
a direct consequence of the fact that Q, Q' have the same 
lower blocks C, D. Now, if the latter of the two possi- 
bilities is the case for some io, the ioth diagonal entries 



of Q{9) and Q'{9) must be different, since Q'{9)igig = 

^M'U = ^M^)^o+^l{o:^'U = Q{o)^.^. + 1, with 

the ioth column of 9'^. As both Q{9) and Q'{9) 
have zero diagonal, this yields a contradiction and we 
have proven the uniqueness of A and B. To prove exis- 
tence, note that for every i, there are exactly two couples 
Oi hi 

above argument that we can always tune {ai,bi) such 
that {A9 + B){C9 + D)~^ has zero diagonal, where we 
take An = ai and Bu = hi for i = 1 , . . . , n. Since 

An Bii 



{a,,bi) s.t. 



each 2 x 2-matrix 



is invertible. It follows from the 



suit. 



A B 
C D 



is invertible, the matrix 
is an element of C', which proves the re- 



□ 



IV. LOCAL CLIFFORD OPERATIONS AS 
GRAPH TRANSFORMATIONS 



In this section, we investigate how the transformations 
(PI can be translated as graph transformations. First we 
need some graph theoretical notions: two vertices i and 
j of a graph G — (V^E) are called adjacent vertices, or 
neighbors, if {i,j} G E. The neighborhood N{i) C V 
of a vertex i is the set of all neigbors of i. A graph 
G' = {V',E') which satisfies V C V and E' <Z E is 
a subgraph of G and one writes G" C G. For a subset 
A CV oi vertices, the induced subgraph G[A] C G is the 
graph with vertex set A and edge set {{i,j} G E\i,j G 
A}. If G has an adjacency matrix 9, its complement C?'^ 
is the graph with adjacency matrix 9 + 1, where I is the 
n X 71-matrix which has all ones, except for the diagonal 
entries which are zero. 

Definition 1: For each i — 1, . . . , n, the graph transfor- 
mation gi sends an n- vertex graph G to the graph gi{G), 
which is obtained by replacing the subgraph G[N{i)], i.e. 
the induced subgraph of the neigborhood of the zth vertex 
of G, by its complement. In terms of adjacency matrices, 
gi maps 6* G 9 to 

g,{9) = 9 + 9A,9 + A, 

where A^ has a 1 on the ith diagonal entry and zeros 
elsewhere and A is a diagonal matrix such as to yield 
zeros on the diagonal of gi{9). 

The transformations gi are obviously their own in- 
verses. Note that in general different gi and gj do not 
commute; however, if 6* G 9 has 9ij = 0, it holds that 
9i9j{^) — 9j9i{^)i one can easily verify. 

Example: Consider the 5-vertex graph G whith adja- 
cency matrix 9ij = 1 for all i ^ j and 9ii = for all 
i (i.e. the complete graph), which is the defining graph 
for the GHZ state. The application of the elementary 
local Clifford operation gi to this graph is shown in Fig. 




FIG. 1: Application of the graph operation gi to the GHZ 
graph. 



1. 

The operations gi can indeed be realized as local Clif- 
ford operations This is stated in theorem 2 and was 
found independently by Hein et al. 5]. 

Theorem 2: Let gi be defined as before and 9 E O. Then 
g,(6) = Q,{6), 

where 



I 

A,; 



diag{9i 
I 



where diag{0i) is the diagonal matrix which has Otj on 
the jth diagonal entry, for j — 1, . . . ,n. 

Proof. The result can be shown straightforwardly by cal- 
culating Q,{e) = {0 + dia,g{e,)){A,9 + I)"^ and noting 
that the matrix Ai6 + 1 is its own inverse for any 6. □ 

The remainder of this section is dedicated to proving 
that the operations gi in fact generate the entire orbit of a 
graph state under local Clifford operations, that is to say, 
two graph states \G), \G') are equivalent under the local 
Clifford group iff there exists a finite sequence , . . . , g.^^^ 
such that gij^ . . . gi^{G) = G' . This result completely 
translates the action of local Clifford operations on graph 
states into a corresponding action on their graphs. In 
order to prove the result, we need the following lemma. 

Lemma 2: Define the matrix class T C Zj^" by 

T = {Ge + D\ G e and G,D are diagonal 
and G9 + D is invertible} 

and consider an element R G T. Choose 9, C, D such 
that R = G6 + D. Define the transformation fi of T 
by fi{X) = X{K,X + X,,Ai + /), for X gT, and denote 
fjki') ■= fjfkfji')- Then (i) there exists a finite sequence 
of fi 's and fjk 's such that 



fjMkM ■ ■ ■ fjikifiN ■ ■ ■ fii {R) — I, 



(4) 



where all the indices in the sequence are different; (ii) 



there exists a unique Qq ~ 
e dom{Qa) and 



C D 



G C', such that 



where gjk{-) := gjgkgji-)- 

Proof: First, straightforward calculation shows that fi 
maps the class of matrices of the form CO + D to itself. 
Furthermore, for each X E T the matrix AiX + XiiAi + 1 
is invertible, which implies that fi maps invertible ma- 
trices to invertible matrices. Therefore each fi is indeed 
a transformation of T. Now, statement (i) is proven by 
applying the algorithm below, where the idea is to suc- 
cessively make each ith row of R equal to the zth can- 
nonical basis vector ef — [0 ... 1 ... 0] , by applying 
the correct /j's in each step. The image of R throughout 
the consecutive steps will be denoted by the same symbol 
R = [rij ) . Now, the algorithm consists of repeatedly per- 
forming one of the two following sequences of operations 
on R: 

Case T. li R has a diagonal entry r^pip = 1 (and the 
«oth row of R is not yet equal to the basis vector e^), 
apply /if,. It is easy to verify that, in this situation, fi^ 
transforms the i^ih row of R into the basis vector e^. 

Case 2: If the conditions for case 1 are not fulfilled, 
apply the following sequence of three operations: firstly, 
fix a jo such that r^,, = and apply fj„ . It can easily be 
seen that then diag(f^) — > d\a,g{R)~\- Rj„^ where R^^^ is the 
joth column of R and d\&g{R) is the diagonal of R. Since 
R is invertible, Rj^ has some nonzero element, say on the 
/coth position r^^j^ = 1. Therefore, the application of fj^ 
has put a 1 on the k^Xh diagonal entry of the resulting 
R. Now we apply /fc(,, turning the k^th row into e^^ as 
in case 1. Furthermore, this second operation has put 
rjgko on the joth diagonal entry and, from the symmetry 
of ^, it holds that rj^ko — fkojo — 1- Therefore, by again 
applying we obtain an ej^ on the joth row of the 
resulting R. Finally, we note that after performing this 
sequence of operations, we end up with an R which will 
again satisfy the conditions for case 2. 

Repetition of these elementary steps will eventually 
yield the identity matrix, which concludes the proof of 
statement (z). 

Statement (ii) is proven by induction on the length 
of the sequence of /i's and fjk's. As it turns out, the 
easiest way to do this is to consider the /i's and fjkS as 
two different types of elementary blocks in the sequence 
The proof will therefore consist of two parts A and 
B, part A dealing with the fi's and part B with the fjk^s. 

Part A: in the basis step of the induction we have 
fi{R) — I, where i? S T. Any such R satisfies R = 



9jMkM ■ ■ ■ gjikiQiN ■■■9ii (^) - Qo{(^), 



AiR + Riiki + / and therefore must be of the form 



R ■ 



Xi+l 

1 



1 



C0 + D 



for some Xj. Then any 9, C, D which satisfy R 
must satisfy Oij = Xj and D = I; moreover 
for j ^ i then the jth row of 9 must be equal to zero and 
if Cii — the ith row of 9 must be to zero. It is now easy 



if = 1 



to see that gi{9) ^ Q[9), with Q = 



C D 



In the induction step, we suppose that the statement 
holds for all sequences fi-^, . . . , fi^ of fixed length N and 
prove that this implies that the statement is true for se- 
quences of length + 1. We start from the given that 

liN ■ ■ ■ fiifi{R) = I 

for some fi, fi-^, . . . , fi,^ and R £ T and we choose 9, C, D 
such that C9 + D = R. Note that it follows from case 
1 in the algorithm in part (i) of the lemma that we may 
take di = 1 = Da, as a single fi (as opposed to an fij) is 
only applied when Ru — 1. Furthermore, we will denote 
by U! the set of all j e {1, . . . , n} such that Cjj = 1. As i? 
is invertible, this implies that D^k = 1 for k G oj'^. Now, 
denoting R' := fi{R), we have /»„... /„ (R') = /, which 
allows us to use the result for length N: for any 9' ,C', D' 
that satisfy C'9' + D' = R' there exists a Q' € which 
has C',D' as its lower blocks such that 9' G dom((3') 
and 

9i.■■■9^A^')^Q'iO')■ (5) 
We make the following choices for 9' ,C' ,D': 

9' = g,{9) 
C = C + Ai 
D' = Z? + diag(6l,)„ 

where diag{9i)^ is the diagonal matrix which has 9ij on 
the jth diagonal entry if j G a; and zeros elsewhere. This 
choice for 9', C , D' indeed yields C'9' + D' = R'; we will 
however omit the calculation since it is straightforward. 
Now, using the definition of 9' and Theorem 2, equation 
(O becomes 



9i 



■ g^,g^{e)^{Q'Q^m■ 



(6) 



It is now again straightforward to show that Q := Q'Qi 
has C and D as its lower blocks. Uniqueness follows 
from lemma 1. This proves the induction step, thereby 
concluding the proof of part A. 



Part B: The proof of this part is analogous to part A, 
though a bit more involved. The basis step now reads 
fjk{R) — I- Now case 2 in the algorithm in part (i) of 
the lemma implies that Rjj = = Rkk and Rjk = 1, as 
only if this is the case, fjk is applied in the algorithm. 
For simplicity, but without losing generality, we take i = 
l,j = 2. Then R must be of the form 



R 



1 

1 







91 
9J 



where the 0i are {n — 2)-dimensional column vectors. 
Choosing 9, C, D s.t. R = C9 + D, the matrix 9 must 
satisfy 



1 

1 



9f 

9^ 



71 02 



where </) is a symmetric {n — 2) x (?i — 2)-matrix with 
zero diagonal; furthermore Dn = = D22, Djj = 1, 
Cii = 1 = C22 and if Cj+2j+2 = 1 then 9ij = 92j = 
(f>kj — 0, for j, fc = 1, . . . , n — 2. We will give the proof for 
Cj+2.j+2 — 0, the other cases are similar. Thus, we have 
to show that there exists a. Qo E C' with lower blocks 
C,D s.t. 912(9) — Qq{9). To prove this, we use theorem 
2, yielding 



Qo 



Ai / 



/ diag(ffi(0)2) 
A2 / 



/ diag(0i) 
Ai / 



where gi{9)2 = (1, 0, 6I13-I-6I23, ■ ■ • ,6'in+02n) is the second 
column of 9i{9). A simple calculation reveals that 



Qo 



A2 









C D 



Ai + A2 / + Ai 

which proves the basis of the induction. 

In the induction step, we again follow an analogous 
reasoning to part A: we suppose that the statement is 
true for sequences fjj^ki , • ■ ■ , /jNfejv of length N and prove 
that this implies the statement for length iV + 1. Our 
starting point is now 

/jjvfew ■ ■ ■ fjikifjk{R) = I 

for some fjk,fnki /jnAn and i? e T. Note that again 
we have Rjj = = Rkk and Rjk = 1 as in the basis step. 
As from this point on the strategy is identical as in part 
A and all calculations are straightforward, we will only 
give a sketch: first we denote R' = fjk{R)', for 9,C,D 



s.t. i? = + D, we define 

e' = 9jk{e) 
C = C + Aj+Afe 
D' = D + Aj+ Ak 

It then straightforward to show that R' = C'O' + D'. The 
induction yields a, Q' £ with lower blocks C", D' such 
that 

93NkN ■ --gjlkiid') = Q'{o'). 

Using theorem 2, we calculate Qjk s.t. gjk{d) = Qjk{S)- 
Then 

93NkN ■ ■ ■ 9nki9jk{0) = {Q'Q]k){0) 

and a last calculation shows that Qq := Q'Qjk has lower 
blocks C and D. Uniqueness again follows from lemma 
1. This proves part B of the lemma. □ 

The main result of this paper is now an immediate 
corollary of lemmas 1 and 2: 

Theorem 3: Let G 6. Then the operations gi, . . . 
generate the orbit of 9 under the action 0) of the local 
Clifford group C' . 



Proof: Let Q 



A B 
C D 



e C' such that 6 e dom((3). 



Now, as C6 + D is an invertible element of T, lemma 2(ii) 



can be applied, yielding a unique Qq 



Ao Bo 
C D 



and sequence of g^'s and g^fe's such that 6 G dom((5o) 
and 



In a final note, we wish to point out that testing 
whether two stabilizer states with generator matrices 
S, S' are equivalent under the local Clifford group is an 
easily implementable algorithm when one uses the binary 
framework. Indeed, one has equivalence iff there exists a 
QeC s.t. 



S^Q^PS' = 0, 



(7) 



as this expression states that the stabilizer subspaces gen- 
erated by the matrices S' and QS are orthogonal to each 
other with respect to the symplectic inner product. Since 
any stabilizer subspace is its own symplectic orthogonal 
complement, the spaces generated by 5" and QS must 
be equal, which implies the existence of an invertible R 
s.t. iS" = QSR. Equation lO is a system of linear 



equations in the 4n entries of Q 



with n ad- 



A B 
C D 

ditional quadratic constraints An Da + BnCu — 1 which 
state that Q £ C''; these equations can be solved numeri- 
cally by first solving the linear equations and disregarding 
the constraints and then searching the solution space for 
a Q which satisfies the constraints. Although we can- 
not exclude that the worst case number of operations is 
exponential in the number of qubits, in the majority of 
cases this algorithm gives a quick response, as for large 
n the system of equations is highly overdetermined and 
therefore generically has a small space of solutions. Note 
that, when equivalence occurs, the algorithm provides an 
explicit Q wich performs the transformation. 



VI. CONCLUSION 



gjMkM ■ ■ ■ 9jiki9iN ■■ - gii {0) - Qoid), 

As Q and Qq have the same lower blocks C and D and 9 
is in both of their domains, it follows from lemma 1 that 
Qo — Q and the result follows. □ 



DISCUSSION 



In this paper, we have translated the action of local 
unitary operations within the Clifford group on graph 
states into transformations of their associated graphs. 
We have shown that there is essentially one elemen- 
tary graph transformation rule, successive application of 
which generates the orbit of any graph state under the 
action of local Clifford operations. This result is a first 
step towards characterizing the LU-equivalence classes of 
stabilizer states. 



The result in Theorem 3 of course facilitates gener- 
ating the equivalence class of a given graph state un- 
der local Clifford operations, as one only needs to suc- 
cessively apply the rule to an initial graph. Note that 
the lemma 2 implies that one only needs to consider se- 
quences gj^kM ■ ■ ■ 9nki9iN • • ■ 5n of hmited length. Fur- 
thermore, the translation of the operations © into se- 
quences of elementary graph operations gets rid of an- 
noying technical domain questions. It is important to 
notice that we have not proven that each Q corre- 
sponds directly to a sequence of g^'s, since, in theorem 3, 
the decomposition into giS depends both on Q as well as 
9. 
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